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Abstract 

The vacuum energy of a scalar field in a spherically symmetric back- 
ground field is considered. Based on previous work [Phys. Rev. D 53 
(1996) 5753], the numerical procedure is refined further and applied to 
several examples. We provide numerical evidence that repulsive potentials 
lead to positive contributions to the vacuum energy. Furthermore, the 
crucial role played by bound-states is clearly seen. 

Pacs number(s): 11.10-z; 02.30.-f; ll.lOGh 



1 Introduction 

The calculation of vacuum energies in the context of quantum field theory under 
external conditions plays an important role in several areas of modern theoret- 
ical physics. External conditions present may be classical fields (for example 
gravitational U and electromagnetic fields, monopoles 0, [|, Sphalerons H or 
electroweak Skyrmions [0]), or they may be realized through boundaries at which 
the quantum field has to satisfy certain boundary conditions (see for example 
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H). External fields most naturally appear as classical solutions of some physical 
nonlinear equations, but they may also be viewed as a localized property of space 
representing a distribution of matter within which the quantum field exists. For 
external fields realized as boundaries this situation is usually called Casimir ef- 
fect. There is a number of generalizations towards weakened boundaries, e.g., 
by introducing semi hard walls |7]] or transparent boundaries realized by delta 
functions f|. For actual interest in such calculations see, e.g., P, |10|. 



A formalism how to deal with arbitrary background fields, which we use here, 



has been developed in |ITJ for a background depending on one cartesian coordi- 
nate and in for a general spherically symmetric background. However, these 
methods had been applied to a square well potential only where the interesting 
quantities can be expressed in terms of Bessel functions. 

The aim of the present work is to show that this formalism is well suited for 
arbitrarily shaped spherically symmetric background fields. We investigate in 
detail different types of background fields as examples. It is observed that the 
main feature of the vacuum energy is determined by the number (and depth) of 
bound states. As long as these are kept fixed the general behavior of its energy 
when varying parameters is the same within the class of examples considered. At 
the same time we find that the presence of bound states does not fix the sign of 
the ground state energy. 

Similar calculations, however with a less extend of details and using other 



techniques, have been performed earlier, e.g., in Jl3| and |L4 . 

The paper is organized as follows. In the next section the basic notations of 
the method used are explained. In the third section the 'asymptotic' contribu- 
tion is calculated. In the fourth section the numerical procedure is developed in 
detail followed by a section containing the examples calculated. The results are 
discussed in the conclusions. 



2 The model, its renormalization and some ba- 
sic results 

We consider the theory described by the Lagrangian 

L = -$(□ — M 2 — A$ 2 )$ + ]-<p(U -m 2 - A'$ 2 )y?. (1) 

Here the field $ is a classical background field. By means of 

V(x) = A'$ 2 (2) 

it defines the potential in ([]]) for the field (p(x), which should be quantized in the 
background of V(x). 
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Our interest is in the vacuum energy of the theory which is defined as half 
the sum of the one particle energies of the field tp(x) 

E,m = lJ2( x h+ m2 y /2 ~^ 2s - ( 3 ) 

(n) 

Here s is the regularization parameter with s — > in the end and \x is an arbitrary 
mass parameter. The \ 2 n ^ are the one-particle energy eigenvalues determined 
through 

(-A + V{x))cj> {n) {x) = Xf n) <f> {n ){x). (4) 

For the moment a finite volume is assumed to have discrete eigenvalues. 

The spectral function appearing in Eq. (H) is related to the zeta function of 
the wave operator (|J), 

Cv00 = 5>? n) +m a )-, (s) 

(n) 

by 

E v m= l -C, v {s-l/2)^. (6) 

The ultraviolet divergences contained in the ground state energy ([3]) needs to be 
renormalized. The divergences are known and best expressed in terms of the heat 
kernel coefficients related to the operator defined in Eq. (H). Then the divergent 
part of the ground state energy reads 

v 1 1 64ir 2 Vs m 2 2/ 

m 2 (1 , V \ „ 

i i + .„^- 2 U, 

with the coefficients Ai = — J dx V(x), A 2 = \ J dx V(x) 2 and A is the spatial 
volume. The contribution of Aq, because it does not depend on the background 
will be dropped from here on without further comment. 
The renormalization procedure is set by the definition 

E™ = E H> [*\-E*°[*\ (8) 

and the corresponding counter terms can be absorbed into a redefinition of the 
parameters M 2 and A of the classical background field whose energy reads 

1 C M 2 A 

£ c iass[<£] = 2 J d£(V$) 2 - —Ax + — 2 A 2 . (9) 
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It might be useful to note that among the divergent contributions (|7|) to the 
ground state energy there are no terms proportional to the kinetic part in the 
Lagrangian of the background Within the given framework this is a simple 

consequence of the structure of the heat kernel coefficients. Sometimes for the 
divergent part of the ground state energy the first 2 terms of the perturbative 
expansion of -Egt^] with respect to powers of the background field $ is taken 
(e.g., in |14]]). This has the disadvantage of containing nonlocal contributions (in 
opposite to the heat kernel expansion expansion which is local) and invites to 
discuss the kinetic term among the divergent ones which is in fact unnecessary. 

As the renormalization procedure is not unique, a normalization condition 
must be imposed. We require 



lim E r v en [<$>] 



0. 



(10) 



For a massive field this is the natural condition requiring the energy of the quan- 



tum fluctuation to vanish for the field becoming heavy. As discussed in [15| this 
delivers a unique definition of the ground state energy (in opposite to the case of 
a massless field, see the discussion in [Hj and also in ||T6f). 



To proceed we follow [|12|] and use the representation of the regularized ground 
state energy in terms of the Jost function fi(k) of the scattering problem corre- 
sponding to (H) 



m 2 } 1 *-* ^fi(ik). (11) 



We remind the definition of the so-called regular solution <p nt i(r) of the radial 



wave equation (see, e.g., fT7[l ) 

1(1 + 1) 



dr 2 



- V{r) + k 2 



<f>n,i(r) = 



(12) 



It is the solution, which for r — > behaves as the free (i.e., without potential) 
solution. The Jost function (and its complex conjugate) are defined as the coef- 
ficients in the asymptotic expansion of this solution as r — > oo: 



fi(p)h l (pr) - f*(p)h + (pr) 



(13) 



where hj(pr) and h~l(pr) are the Riccati-Hankel functions. 

Instead of dealing explicitly with eigenvalues, the vacuum energy is expressed 
in terms of the Jost function. Although representation ( |TT| ) does not reveal explic- 
itly the dependence of the energy on the bound states, in |L2| it was shown that 
their contribution is correctly included. This, and the non oscillating behavior of 
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the Jost function on the imaginary axis (cf. (|ITD) gives this representation cal- 
culational advantages with respect to alternative representation (e.g. integrating 
the scattering phases over the real fc-axis). 

To proceed we have to perform the analytic continuation in s to s = in 



representation 



of E™ n using 



PI) for i^[<I>]. Although we known that the 



result is finite, this task requires some work. The point is that one cannot put 
s = under the sign of the sum and the integral in 



Tl|) being inserted into Eq. 



(||p. The reason is the lack of convergence for k and I large. Therefor we use the 
uniform asymptotic expansion In f^ ym (ik) of the logarithm of the Jost function 
for both, k and / large. We subtract and add it to the integrand. As a result we 
get the sum of two parts 



EFffl=E f [*) + E„[$] 



with 



oo 00 

71 1=0 L 



(14) 



ik)\ (15) 



and 



COS! TCS 



M^^(i + i/2) I dk [k 2 
1=0 £ 



m^-^lnfr m m- < v .(16) 



The first contribution, is called the finite part. There the regularization 

parameter s has been put to zero directly because the sum and the integral are 
convergent due to the subtraction. The second contribution, E as [$] , is called the 
asymptotic contribution. It has, as we know, a finite limit for s —>■ which will be 
determined below. The asymptotic expansion In fi Sym {ik) for v — > 00, k — > 00, 



vjk fixed, was determined in JT2[ up to the required order of v 3 . It reads 

00 

r V(r) 



\nff sym {ik) = 1- j dr 







1/2 



+ 



16i/ 3 
1 



dr 



3/2 



dr 



r 3 V 2 (r) 

ol 3/2 : 

+ (?) 2 



i + (v) 2 i + 



(17) 



with = I + 1/2. 
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3 The asymptotic contribution 



First of all the /^-integrals in (|T6| ) can be done explicitly using 

oo 

jdk [k 2 -m^^- k 

m 



1 + 



r(n/2) 



This naturally leads to functions of the type 



\mr J 



n-l 

2 



18) 



f(s;a;b)= V) i/» ( 1 + (— Y 

^=1/2,3/2,... v 



-s-6 



(19) 



which we need around s = —1/2. The way to deal with this functions has been 
developed in [Tj| and for completeness the relevant results are provided in the 
Appendix, (PI)). Using these, it is a simple exercise to find 



Kt n m 



57T 



dr r V (r) ln(mr) 



(20) 



2 poo roo 

— / dr V(r) / dv — 
2vr Jo K) J 1 



V 



8vr Jo 



dr 



+ e 



rV 2 (r) - -\^(r) 



2nu 



In |z/ 2 — (mr) 2 | 



:x3 / d 
dp , 

di/ 1 + e 2 ™ 



In v 2 



[mr] 



dr V(r) / dv 
87T ./o Jo 



d /Id v 2 



48vr 



o?r V(r) 



dv 



dv \vdvl + e 2nu 
d ( \ d \ d v A 



dv \ v dv v dv 1 + e 27r 



In \ v 2 — (mr) 2 | 
In I v 2 - ( 



mr 



suitable for numerical evaluation. 



4 Numerical determination of the Jost function 
and the integration procedure 

Under the assumption that the potential has compact support, let's say V(r) = 
for r > R, the regular solution may be written as 



<M0 = H P (r)Q(R - r) + - [fi( P )K(pr) - / { *(p)A+(pr)j 0(r - fl). (21) 
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This solution and its derivative must be continuous at r = R. So the matching 
conditions are 

ui, p (R) = l[fi(p)K(pR)-ft(p)ht(pR) 
u[ p (R) = l - P [fi(p)K'{pR)- K(p)ht'(pR) 

From these conditions we get the Jost function in terms of the regular solutions, 
flip) = ~ (pu^h+'ipR) - u[ p (R)h+(pR)) , (22) 

where we used that the Wronskian determinant of hf is 2i. 

The assumption of the potential V(r) having a finite support is natural within 
the framework of a numerical calculation. Anyway the potential is assumed to 
vanish for r — ► oo. So, after some value of the radius it can be viewed as zero 
within the required precision. This simplifies the numerical procedure consider- 
ably. Another situation would occur for a Coulomb like potential which we do 
not consider in this paper. 

In order to determine a unique solution of the differential equation fll~2]) , we 
need to pose the initial value problem corresponding to 4>i, p (r) being the regular 
solution. We have for r — > 

Mr) = u ltP (r) ~ Uw) ~ r{l f 3/2) (f^T ■ 
We change to the new function 



i \ v /7r fpr\ l+1 . 

with the initial value gi, p (0) = 1- The differential equation for gi tP (r) reads 

or, going to the imaginary p-axis and writing it as a first order differential equation 
with (d/dr)gi !ip (r) = v hip {r), 



± ( gi,i P (r) \ _ ( 1 \ ( gi jip {r) 

dr I v hip (r) i I V(r) + p 2 -Ml + 1) M v ljip (r) 



(23) 



To fix the solution uniquely we only need to fix f; i i P (0). A power series ansatz 
for gi/ ip {r) about r = shows that for V(r) = 0(r~ 1+e ) the condition reads 
v i,ip(fy — 0- With this unique solution of (p3|), the Jost function takes the form 



fiiip) 



T(/ + 3/2) V 2 



g ljip (R)K l+3/2 (pR) + l -gl ip (R)K l+1/2 {pR)} .(24) 



Finally, doing a partial integration and the substitution q = \Jk 2 — m 2 in Eq. 
(|16"D, the starting point for the numerical evaluation used for Ef is 



1=0 



I — r 

= + 1/2) / <k ^fi{iy/¥+ 

7T In L 



in/; 



asym 



(25) 



where In J z asym is given by Eq. (|ITD . 



5 Examples 

The method described above works for radially symmetric potentials $(r) having 
compact support $(r > R) — 0. In order to have a finite classical energy (Eq. |9|) 
we also demand <3>'(r = 0) = $'(r = R) = 0. 

In the following we use dimensionless quantities 

(26) 



e = 


ER, 


^ = 


mR, 




r 


P = 


r: 


V(r) = 


A'$ 2 (r) 











such that e = e[(p, ft]. Therefore, large /x corresponds to large i? and to large m, 
respectively. 

We studied two types of potentials. Type A is lump-like concentrated around 
r = whereas type B has its maximum at R/2, acting like a spherical wall: 

, i a(l — p 2 ) 2 , , 

Mp) = . 2 » ( 2? ) 

a + p 2 
16ap 2 (l-p) 2 

<Pb(p) 



a + (1 - 2p) 2 ' 

The parameter a allows to vary the shape of the potentials, as indicated in the 
insets of the figures. 

Our numerical implementation of Eq. ( p5|) showed that the sum over angular 
momenta converges quite fast, usually after 10-50 terms. Hereby, the numerical 
g-integration procedure has to solve the system of ordinary differential equations 
Eq. fl2~3D for every evaluation of the Jost function fi(ip). 

As a check, we compared our results with the asymptotic formulae 



e ~ 




A 2 log - for u — >• . , 

& ^ P (28) 
for p — > oo , 
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where A2 = \ j dxV(r) 2 is the second heat kernel coefficient and 
A3 = I f dx ^V(r)AV(r) — V^(r) 3 )) is the third one. These relations follow 
from scaling arguments together with the well known property of the heat kernel 
expansion being the asymptotic expansion for large mass. 

The most determining property of the potential V{r) in the attractive case 
A' < is the number and depth of bound states. The parameter dependence of 
the I = bound states, determined as zeroes of the I = Jost function, Eq. (|24]) , 
is shown in figs. [l| and [| Number as well as depth of the bound states increases 
with increasing |A'| and with increasing a, too. 

As Eq. ([25|) shows, the Casimir energy develops an imaginary part (corre- 
sponding to particle creation by the external field) when the g-integration passes 
a zero of the Jost function. Ef is real as long as fi is larger than the depth of 
the highest bound state. Therefore we indicate in our figures |3|-^]the positions of 
bound states on the fi axis by vertical lines where the plots of the Casimir energy 
end. 

To show clearly the influence of the form and spectrum of V on the vacuum 
energy, we consider the dependence of the energy for fixed A' and different values 
of the parameter a, figs. ^| and |], as well as for fixed shape a and different depths A' 
in figs. |5] and |[ 

These figures demonstrate that the vacuum energy is determined by two com- 
peting contributions, the positive contribution of the scattering states and the 
negative contribution of the bound states of the potential. 

Let us start with a description and interpretation of figs. |^ and [|. For a large, 
the contribution of the bound state (s) to the vacuum energy is large enough 
to overcompensate the contribution of the scattering states and the energy is 
negative. At some "critical value" of a, the energy of the only remaining bound 
state is so small that the positive contributions of the scattering states get the 
upper hand, and the vacuum energy is positive. Finally, for a —>■ we have 
normalized _E£, en [$ = 0] = 0, such that at some point E^ n starts to decrease 
again with decreasing a. 

The same features are clearly recovered in figs. |^ and |6|. For large enough 
(—A') a negative vacuum energy energy is obtained. With decreasing (—A') the 
energy increases and can be positive or negative depending on the parameter \i. 
Again, at some point the scattering states dominate, the energy gets positive and 
tends to zero as A' tends to zero. For positive A' the energy starts to increase 
again, being virtually identical for the two values A' = ±2, and it seems that E^ en 
is a monotonically increasing function of A' for A' > 0. Since we see no qualitative 
difference between our two potential types A and B, we expect these features to 
be universal. 

Given that the dependence of the energy on A' is somewhat nontrivial, let us 
show this dependence explicitely for type A potentials in fig. |7|. 

As long as no bound states exist, the sign of A' is not important and the 
energy is nearly symmetric about A' = 0. The contribution to the energy of 
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the scattering states does not depend much on the sign of A'. The symmetry 
is disturbed as soon as a bound state shows up and the energy then decreases 
with decreasing A' because the bound state depth as well as the number of bound 
states increases. 



Conclusions 

We demonstrated that our method is quite feasible for calculations of the Casimir 
energy in realistic 3-dimensional examples for arbitrary spherical symmetric ex- 
ternal potentials as long as they are smooth and have compact support. 

The behavior of the Casimir energy shows some interesting features. The 
existence of bound states is a necessary but by no means sufficient condition for 
a negative vacuum energy. Even in the presence of bound states the sign of the 
energy can depend on the scale parameter /i, i.e the mass of the quantum field 
measured in units of the size R of the potential. 

For shallow potentials without bound states (small |A'|) the Casimir energy 
is nearly invariant under a change of the sign of the potential whereas for deep 
attractive potentials the negative contribution of the bound states dominates the 
vacuum energy. 

On the other hand we found no qualitative differences between a lump-like 
and a ring-like classical background configuration. 

The method developed here can be used for the calculation of quantum cor- 
rections to classical objects like Skyrmions and others. 

KK had been supported by DFG under grant number Bo 1112/4-2. 



Appendix: Analytical continuation of the asymp- 
totic contributions 



The essential formulas for the basic series f(s;a;b), Eq. ([19]), in the calculation 
are the following (we use x = mr): 

^=1/2,3/2,... V 





oo 



v 2n+l / 2 \~s 



■ 
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i/=l/2,3/2,. 



z/ 2n | 1 + ( - ) 



1 r(n + l/2)r(a - n - 1/2) , n 



-l) n 2sin(7rs) / dv 



Vis) 



,2» 



+i 



(30) 



X/ 



2 

) -1 



1 + e 27W 

Using partial integrations one can get representations valid for all values of s 
needed. The explicit results applied for the calculation of ( |16D are [we shall use 
the notation f(a; b) = /(— 1/2; a; £>)], 



ds 



/(i; i/2) 

W(s;l;l/2) 



/(i; 3/2) 



1 2 1 

— H 

2 24 



1 



-x 



x 



dv- 



1 + e 



2tt!/ 



hi 



V\ 2 



2(s + l/2) 



+x z 



du 



X 111! 



d 1 



/(3;5/2) 



2(s + 1/2) 



+■ 



4 />oo 



/(5;7/2) 



2 Jo 

x e 



dv 



dvl + e 2nv 
+ (In a; - l/2)x 4 
d Ad v 2 



In z/ — 



2(s + 1/2) 
+ — / dv 



dv V v dv 1 + e 2?r!/ 
+ (lnx - 3/4)x 6 

d f 1 d 1 d v A 



In i/ — x 



<iz/ \v dv v dvl + e 2nu 
which provide the needed expansions about s = in (IB). 



(31) 

In \v 2 — x 2 \ , 
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Figure 1: Energy of bound states for type A potential and negative A'. (The inset 
shows ip(p) where V(p) = \'(p 2 (p)) 
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Figure 2: Energy of bound states for type B potential and negative A'. (The inset 
shows (p(p) where V(p) = \'(p 2 (p)) 
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Figure 4: Energy for type B potentials of different shapes with A' = —100. 
positions of bound states (BS) at the ji axis are shown as vertical lines. 
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Figure 7: Energy for type A potentials for different magnitudes A' and equal 
shape parameter a = 1 as a function of the parameter /i. 
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